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Charge fluctuations in nano-circuits with capacitor components are shown to give rise to a novel 
type of long-ranged interaction, which co-exist with the regular Casimir/van der Waals force. The 
developed theory distinguishes between thermal and quantum mechanical effects, and it is applied 
to capacitors involving graphene nanostructures. The charge fluctuations mechanism is captured 
via the capacitance of the system with geometrical and quantum mechanical components. The 
dependence on the distance separation, temperature, size, and response properties of the system 
shows that this type of force can have a comparable and even dominant effect to the Casimir 
interaction. Our results strongly indicate that fluctuations induced interactions due to various 
thermodynamic quantities can have important thermal and quantum mechanical contributions at 
the micro- and nanoscale. 
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I. INTRODUCTION 


Fluctuation-induced interactions have widespread applications in materials and micro- and nano-scaled devices. 
The much studied Casimir and van der Waals interactions are due to electromagnetic modes fluctuations captured via 
the dielectric and magnetic response properties of the objects[Tl [5]. Such forces are universal, and they are especially 
prominent at micron distances and below. Fluctuations of many other observables are also possible, which may give 
rise to different interactions [3]. In particular, voltage fluctuations in capacitor systems |4] and wires[5] have been of 
much interest, especially for the operation of devices |^. Charge fluctuations and the induced forces are of relevance 
to biological and chemical matter. Specifically, it has been shown that thermal charge fluctuations in ionic solutions 
can generate an attractive long-ranged dispersive force even between molecules charged with the same sign 0,El- 
Fluctuations originating from charge disorder on neutral slabs have been shown to give an additional contribution to 
the net interacting force while completely masking the typical Casimir-van der Waals interaction!^. 

The isolation of single graphene layers and synthesis of related nanostructures, such as carbon nanotubes and 
graphene nanoribbons (GNRs), have brought new directions in electromagnetic fluctuation induced phenomena. For 
example, the Casimir/van der Waals force involving graphene systems has non-trivial scaling behavior, which can be 
tuned via temperature, chemical potential and doping modifications [101112) . Understanding the electromagnetic fluc¬ 
tuations in the context of building a fundamental understanding and making technological designs for such materials 
can hardly be overestimated. 

Capacitor systems involving graphitic nanostructures have also been studied extensively. The application of a bias 
on a graphene system above a substrate is of great relevance to many quantum devices. A key factor has been the 
control of the local graphene electrochemical potential for the device functionality. Capacitance measurements can 
also be utilized as means to probe the basic electronic properties of graphitic nanostructures and to learn about their 
Dirac-like nature. In particular, the 2D character of graphene is exhibited in the quantum capacitance [13],[M], a 
concept describing the penetration of the electrostatic field due to partial screening m- Despite the presence of 
charges and voltage bias in graphene based capacitive systems, their fluctuations and subsequent induced interaction 
effects have never been considered. In this paper, we show that fluctuating charges in a capacitor system are governed 
by fluctuation-dissipation relations, giving rise to a novel long-ranged interaction. The geometrical size of the capacitor 
is found to play an important role, as it can be used as an effective way to modulate this interaction. 

The charge induced dispersive interactions are presented via a general theory utilizing the capacitance concept. We 
distinguish between thermal and quantum mechanical effects as the characteristic dependences on distance, temper¬ 
ature, and other factors are shown. This Casimir-like force is examined in graphene based capacitors and compared 
to the typical Casimir interaction. 



FIG. 1. (Color online) (a) Ideal system under consideration: A parallel plate capacitor connected by a wire supporting voltage 
fluctuations 5V. One of the plates can be a graphene ribbon with a width w above a metal substrate separated by a distance 2 . 
(b) Experimentally realistic setting: The wire is represented by a metallic connection between the “bottom” metal plate and 
the “top” graphene nanoribbon. 


II. CHARGE FLUCTUATION INTERACTIONS AND CAPACITANCE 

The capacitor system under consideration is shown in Fig. 1. It consists of two plates connected by a wire. The 
voltage fluctuations of the wire induce excess fluctuating charges on the capacitor. As a result, a force originating 
from these type of fluctuations is induced between the plates. The description of the interaction is closely related 
to the surface charge density response 6a due to the external voltage SV via the interaction Hamiltonian SH = 
-fd^r'a(r',t}6V(r',t): 
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Sa{r,t) = f dt [ d^r'C{r,r'-t — t')SV{r',t'), (1) 

J —OO J 

A key quantity here is the capacitance of the system, which can be written in terms of the commutator of the 
charge density-density correlations as|16j 

C(r,r';t-t') = ^([<7(r,t),a(r',0]), (2) 

where (...) = Tr(p...) with p being the unperturbed density matrix. 

The density correlation can further be represented as 

]^{5a{Y,t)5a{Y' ,t') + 5a{v' = J ^5(r, r'; (3) 

and the fluctuation-dissipation theorem is utilized via the structure factor 5(r,r';w) = hQmC{r,r']uj)coth ^ t ) • 
Thus the charge fluctuations are given in terms of Matsubara frequencies a;„ = 2TmkBT/h 


{Sa‘^) = kBT ^ lim C(r,r';i|w„|) =/cbT ^ lim / —^C'(k;z|a;„|) 

^' r-i-r' ' r—i-r' / (z7rl 

-JO— r>r~\ rt — r>r~\ t/ \ / 


o*k-(r-r') 


( 4 ) 


{6a^) ^ (Sq^) = kBT ^ C(f|cc„|), (5) 

n— — (x> 

where k is the 2D wave vector. When the voltage fluctuations are distributed uniformly on the plates, the surface 
charge density fluctuations can be written in terms of the total charge fluctuations according to Eq.([^. 

Using the electrostatic energy stored in the capacitor U = q^ j^C , one finds the force per unit area between the 
capacitor plates as 


/ = 


1 9C7-1 
2A dz 


{5q^) 


kBTdC-^ 


OO 

n——oo 
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where A is the area of the plate, and C{(jj) is the capacitance of the system that includes the frequency dependent 
response of the wire. One should note that the above expression does not take into account the frequency dependent 
response of the graphene sheet, for example. Eq.(^ constitutes a general expression for a charge-induced fluctuation 
interaction. This Casimir-like effect is determined by the capacitance of the system. To understand the phenomenon 
further, the theory is applied to graphene systems with emphasis on thermal and quantum mechanical contributions. 

The capacitance is obtained by considering the charging of the capacitor, which leads to redistribution of carriers 
in space in order to minimize the electrostatic energy. Consequently, for classical systems it is determined by the 
geometry of the system via the geometrical capacitance Cq. For thin films, such as graphene and GNRs, the surface 
charges cannot completely shield the electrostatic field. This results in raising of the chemical potential to account 
for the increased density of states necessary for complete shielding m- The effect is associated with the quantum 
capacitance Cq, and the total capacitance C is determined by an addition of capacitors in series according to 


1 _ 1 1 
C ^ A^' 


( 7 ) 


For the rectangular system in Fig. 1 the geometrical capacitance can be found[TS] as 


Cn = - 


2z arctan (-;—) -|-^ln<l-|- 


(£) 


( 8 ) 


One notes that when the width of the ribbon w is large, the parallel plate capacitance is recovered Cq = A/id-nz). 

The quantum capacitance is associated with the particular material. For graphene it may be obtained via the 
random phase approximation m, where the particle density is determined from the density matrix p and the particle 
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density operator - n(r) = Tr [(5/95(r — ri)]. Using the graphene band structure described with a linear two band model 
near the characteristic K point, Uk.s = svohk with vq = 10® m/s, it is found that 

= _2e^ (fk+qs' -/ks)|(5'k + q|ks)p 

A ^ Uk+q.s' — Ek,s — Eo 

k.ss' ’ 



where fi is the chemical potential, s, s' = ±, /k,s = ^-f!)/kgT and (s'k + q|ks) is the overlap integral between 

the wave functions [20]. In the static limit (a;,k) —0, which is of relevance here, one finds [211 [22] 




ge^ksT 


( 2 ^) 


2d{ne - Uh) 
® dg 
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where rig, nu are the electron and hole concentrations respectively. The above expression reflects the fact that there 
is only a hnite amount of charge available on the graphene sheet, Q = e(ne — Uh), due to the reduced dimensionality 
of graphene and its subsequent restriction on the density of states. Using Eqs. an effective penetration distance 

may be defined for graphene as d* = {hvo)^/ [4:ge‘^kBT\n{2cosh{g./2kBT))]. At room temperature with no chemical 
potential, d* « 1 nm. 

The quantum capacitance for GNRs can be obtained in a similar manner. For thin ribbons, however, the energy 
band structure is quantized due to the vanishing wave functions at the edges [20] • The effect of the energy quantization 
in the quantum capacitance can be seen in Fig. 2a for an armchair GNR with re « 12 nm. While the graphene Cq 
monotonically increases as a function of g,, the ribbon Cq experiences peaks |18j . This peak structure is much reduced 
if T is raised to room temperature. 





FIG. 2. (Color online) (a) Cq as a function of chemical potential for graphene and a GNR with w = 12.4 nm. (b) Charge 
fluctuations thermal force normalized by /o = —kBT/{2Az). The area of all ribbon structures is 1 gm^ and g = 0. The 
graphene/metal Casimir force normalized by fo at room temperature is also shown. The metal plasma frequency is taken as 
Ldp = 9 eV. (c) The normalized charge fluctuations thermal force vs chemical potential for different GNRs at T = 300 K and 
separation 0.1 gm. 


III. RESULTS AND DISCUSSION 

The charge fluctuations force is determined by the total capacitance, which contains the frequency dependent 
quantum capacitance according to Fq. §. We distinguish between thermal and quantum mechanical contribution 
in a similar manner as in the case of the typical Gasimir force. The thermal n = 0 Matsubara frequency is considered 
first. The force per unit area for the capacitor with a thicker GNR, in which the discrete band structure is not 
apparent, (Fig. 1) is found analytically 
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For a graphene parallel plate capacitor, the geometrical capacitance is replaced by Cq = A/Airz and w —>■ 00 in /t. 

The characteristic behavior of the thermal charge fluctuations interactions can now be determined using Fq.(ll|. 
Fig.iD shows how the thermal force depends on the separation for both graphene and GNRs with different widths. 
For separations 2 ; >> d*, one finds that Cq » Cq, and the thermal force is determined mainly by the geometrical 
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capacitance. In the case of graphene, the force reduces to /o = — For separations z « d*, Jt is determined 

mainly by the quantum capacitance. For thinner ribbons, the magnitude of the force is reduced, as seen in Fig. 
While the effect of Cq for graphene is relatively small over a large part of the distance range z, its role for thinner 
ribbons can be much enhanced (Fig. 2c). This feature, which reduces the magnitude of /t, is attributed to the 
quantized band structure of the thin ribbons. 

It is important to compare how the typical Casimir force differs from the charge fluctuations induced interaction. 
The Casimir graphene/metallic substrate interaction has been calculated by the Lifshitz approach [TU]. In Fig. 2b we 
show results for an infinite graphene sheet, described by the Dirac model, and a typical metal with a plasma frequency 
of Wp = 9 eV. Due to its 1/z dependence, however, fr is more important for larger z. It should be noted that the 
region where the size of the Casimir force is comparable to the charge fluctuation force is in that separation length 
where the thermal Casimir force is dominant. This reduction in the distance scale where thermal effects take place, 
as compared to other systems, is due[53] to the reduced dimensionality of graphene. 

Fig.H shows how the thermal charge fluctuations force evolves as a function of /i for GNRs with different widths. 
Due to the quantized electronic structure, fi affects the force much more for thinner ribbons. For thicker ribbons, 
however, the interaction is hardly changed upon fi. 

The quantum mechanical contribution to the charge fluctuations force is closely related to the frequency dependent 
total capacitance.One must account for the particular mechanism for charging the capacitor, which happens via the 
connecting wire, as shown in Fig. 1. Taking a Drude model for the wire conductivity leads to a frequency dependent 
total capacitance according to 


1 


1 1 


iujR — 


( 12 ) 


where R is the resistance of the wire and r is the relaxation time mm- Using Eqs. (HU]), the charge fluctuations force 
can now be calculated by taking into account all Matsubara frequencies. Results for the distance and R dependence 
of the force are shown in Fig. It is found that the role of R is more pronounced for smaller z as the effect is stronger 
for smaller R values, according to Fig. §1. Similar conclusion is obtained by considering the dependence of / vs R, 
shown in Fig. If. which indicates that larger resistance (for given t/Cq) leads to thermal fluctuations dominating 
the interaction. 



FIG. 3. (Color online) (a) Charge fluctuations force vs. distance between ribbon of width 0.1 fim and area 1 and a metal 
plate at room temperature, where the C^{u}) is used. Typical values for r and R are taken, (b) The normalized by fo charge 
fluctuations force as a function of the wire resistance for different relaxation times and a plate separation of 0.1 fj,m. 


The quantum mechanical effects related to the mechanism of charging the capacitor can be analyzed further by 
considering the situation when RCq is small compared to r and Cq can be neglected. The total capacitance in this 
case becomes C*{oj) « iSrCot ■ From Eqs. (j^, the force is found as 

/ =- ^ coth ( - ^ ^ . (13) 

AAz^/WJ^ \2kBT^/RC^J ^ ’ 

In the classical limit of large T, one recovers /o = — quantum mechanical limit, where hj^RC ^^t ^ ksT, 

we obtain 


/ = - 


h 


2A^/z Y ARt 


(14) 


Eqs. (13 14) are instrumental in better understanding the / dependence upon z in Eig. 3. As the distance is 
increased, the thermal effects become more important according to the ijz asymptotic behavior. However, for 
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smaller separations, the quantum mechanical contributions become more relevant according to the \pz asymptotic 
behavior. In all cases, these effects are more pronounced for smaller i?, otherwise, the charge induced interaction is 
determined by the thermal fluctuations. 

Casimir and van der Waals forces are universal and of great importance in nature, especially when the scale goes 
down to micrometer distances and below. This type of interaction, induced by dipolar fluctuations, can be essential 
in determining the behavior of nanostructures and devices. Yet there are many other fluctuation induced forces that 
become important at such small scales. Here we consider charge induced fluctuation forces in capacitive systems. The 
developed theory utilizes static and frequency dependent capacitance, and it is applied in graphene based structures. 
The total capacitance contains both a geometrical and a collective electronic component. Just like in the Casimir/van 
der Waals forces, we distinguish between thermal and quantum mechanical effects. The thermal fluctuations are 
pronounced at larger separations where the geometrical capacitance dominates, while the quantum fluctuations are 
relevant for smaller separations. 

The origin of this induced interaction comes from fluctuating charges, which is a key difference when compared to the 
Casimir/van der Waals interaction that originates from fluctuating polarizations. This is important for understanding 
the 1/z asymptotic distance dependence as compared to the dependence of the thermal Casimir force. Similarly, 
the quantum mechanical regime of is of much longer range as compared to the metal/graphene quantum 

mechanical asymptotics m, m- It is also important to note that the geometrical capacitance takes into account 
the geometry of the system in a rather straight forward manner. This is unlike the Casimir force, where non-trivial 
geometry is usually difficult to describe theoretically. This is reflected in the fact that capacitance can be found for 
objects with non-trivial boundary conditions via electrostatic methods, while electromagnetic boundary conditions 
are not easily solved beyond highly symmetric extensions. 

Another interesting feature is the role of the size of the system. The charge fluctuations force is explicitly dependent 
on the area of the interacting objects, which is not the case for the Casimir/vdW force.The charge fluctuations 
interaction is ^ 1/A for the thermal and ^ 1/A^/^ for the quantum mechanical regimes. Balancing the longer ranged 
asymptotics and the area dependence can be used to enhance the charge fluctuations force over the typical Casimir 
interaction. Systems, such as graphene ribbons, which have smaller area as compared to large graphene flakes, are 
especially good candidates to observe this effect. 

In summary, we argue that charge fluctuation forces are always present in capacitor-like systems. They must be 
considered in conjunction with the standard Casimir interaction as the charge fluctuation force can be comparable 
or even bigger in magnitude. Charge fluctuation forces enable further probing of thermal and quantum mechanical 
effects due to fluctuation induced phenomena in nanostructured materials such as graphene and GNRs. We finally 
point out that it is worthwhile to pursue further the effects of fluctuation forces due to not just electromagnetic 
fluctuations but also due to other thermodynamic parameters, for they will necessarily have an important effect in 
many nanostructures and devices. 
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